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Abstract. 

We present mid-infrared (4-18 fim) observations of the microquasar GRS 1915-1-105 obtained with ISOCAM, the 
camera on board the Infrared Space Observatory (ISO), in 1996 April and 1997 October. The first observation 
probably occurred during a flaring event with oscillating synchrotron emission. The 1997 observation occurred a 
few days before a major relativistic ejection, during a plateau state of inverted-spectrum radio emission and hard 
quasi-stable X-ray emission. The K-M giant donor star in GRS 1915-f 105 cannot account for the mid-IR emission 
and we discuss the possible additional components depending on two absorption laws. Thermal emission from dust 
seems unlikely. The flat mid-IR spectrum obtained during the plateau state is likely to be synchrotron emission. 
It would be the first evidence of the infrared extension of the radio synchrotron emission from the compact jets, 
although optically thin free-free emission from an X-ray driven-wind from the accretion disc cannot be excluded. 

Key words. Stars: individual: GRS 1915-1-105 - X-rays: binaries - Infrared: stars - Stars: circumstellar matter - 
ISM: jets and outflows - Stars: winds, outflows 



1. Introduction 

GRS 1915-1-105 is a transient hard X-ray source discov- 
ered i n 1992 with Granat/ WATCH l|Castro-Tirado et alJ 
Il992l [ 1994|) . It is one of the brightest X-ray source 
in the Galaxy, highly variable in intensi ty and activity 
(jHarmon et al.l Il993 iGreiner et alJll996l). Intense flares 
are also observed in radio ijRodrfguez fc Mirabell [1993) 
and led t o the discovery of superlum inal ejections with 
the VLA l)Mirabel fc Rodriguezll99^ i.e. structures mov- 
ing with an apparent velocity higher than the light speed 
on the sky-plane. This first superluminal movement ob- 
served in our Galaxy, very similar to the ones observed 
m quasars like 3C279 ijCotton et al.l ll979V confirmed the 
microquasar nature of GRS 1915-1-105. The microquasar 
term refers to X-ray binaries showing radio (persistent 
or not) jets imitating at much sma ller scales the radio 
lobes of quasars and radio galaxies I Mirabel et al.lll993 

iMirabel fc Rodrfguejll999|) . 

GRS 191 5-H05 has a radi o llMirabel et alJll998b^ and 
an IR (Mirabel et al.lll993a|) count erpart, and is highl y 
variable at all these wavelengths ijMirabel et alJ Il99^ . 
Because of the high interstellar extinction along the line 
of sight to GRS 1915-f 105, the type of the companion 
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star and thus the masses of the bi nary compone n ts were 
difficult to identify. According to iGreiner et all ll2001hf) 
the donor star is of K-M III type correspon ding to a 
late- typ e low mass (~ 1-1.5 M©) giant star. iGreiner et alJ 
l|2001aj) found the following binary parameters: orbital pe- 
riod Porb = 33.5 days, mass function f{M) = 9.5±3.0 Mq 
which is the lower limit to the mass of the compact ob- 
ject, using Afd = 1.2 ± 0.2 Mq for the mass of the donor 
star and 6 — 70° ± 2° for the inclination angle implies 
Mc = 14±4 Mq for the mass of the compact object which 
is the most massive stellar black hole candidate. 



1.1. Multi-wavelength variability 

GRS 19154-105 is unique for its X-ray variability at 
all ti me scales from less than one second to several 
days ijGreiner et alJ Il996(l . RXTE observations of the 
sou rce were cla ssified in different classes (12 classes found 
by iBello ni et al.. 2000a and a new class discovered by 
iKlein-Wolt eraLll2002l) according to their flux and spec- 
tral properties and corresponding to transitions between 
3 basic states. Flaring events seen in radio, near-IR (K- 
band) and X-ray with vario us periodicity and inten sity 
were classified in 3 classes bv lEikenberrv et al.l l(200nl) as 
in the following. 
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The first class is composed of the giant flares that 
were observed in radio and X-ray llMirabel fc Rodrfgued 
1994; Rodriguez fc Mirabell Il999t iFender et all Il999at 
Dhawan et al.. .200Q) and correspond to major ejection 



events followed as super luminal motions in radio images. 
These flares are characterized by very rapid rise time (less 
than 1 day) with an increase in the radio flux density 
by 1 or 2 orders of magnitude (reaching 600-1000 mJy) 
and optically thin e xponential decay lasting ten or so days 
([Foster et al.ll99^ . Continuous short period radio oscilla- 
tions (between optically t hin and thick states) can h appen 
during the flux decrease llFender et alJll999ai Eool . 

These 20 to 40 minutes quasi-periodic oscillations 
form the second type of flares, first ob served i n ra- 
dio w ith the Ryle Telescope at 15 GHz l)Poolevl Il995l 
Il996l) and with the Very L arge Array (VLA) at 3.6 cm 
l|Rodrlguez fc Mirabell 1 1997|l . These oscillations are very 
co mmon as shown by the Ryle Telescope monitoring 
of IPoolev fc FendeJ l|l997l) who associated them with 
soft X-ray dips on ti me-scale of ~ 30 minu tes. The lat- 
ter were discovered bv ICreiner et al.l 



with RXTE 

PC A obs ervations and were interpreted bv iBelloni et al. 
lll997albl) as the repeated disappearance and refilling 
of the inner accr etion disc. N e ar-IR oscillations were 
also observed by iFender et alJ (^£3) who suggested 
that they were the high frequency tail of radio syn- 
chrotron emi ssion. Simultaneous near -IR and X-ray ob- 
servations bv lEikenberrv et al.l l) 199881) revealed the cor- 
respondence between these X-ray dips/flares and near-IR 
(2.2 /xm K-band) flares consistent with synchrotron emis- 
sion from ejected plasma bubbles or "baby jets" under- 
going adiabatic expansion losses. The synchrotron nature 
of these o scillations was c onfirmed by simultaneous near- 
IR/r adio JPender fc Poolev ,.199&) and ncar-lR/radio /X- 
ray l)Mirabel et al.lll998l) observations of these ~ 30 min- 
utes dip/fiare cycles, with wavelength-depen dent time de- 
lays (r adio-radio or radio-infrared). R.ecentlv lFender et alJ 
(120021) demonstrated that this time delay is variable be- 
tween epochs, with a possible correlation with the oscilla- 
tion amplitude, but they still cannot distinguish between 
models of dis crete ejections (as t he Ivan der LaanI Il966l 
model used bv lMirabel et al.(ll998l) or shock s propagating 
along quasi-contin uous flows (jBlandford fc Konigl .1979; 
iKaiser et al]l2000|) . 

The amplitudes of the quasi-periodic oscillations are 
about 40-50 mJy in radio, or ~ 70-80 mJy for the larges t 
and -300 mJy in millimetre (|Fender fc PoolevI \2Q0(h . 
In K-band the obs erved (not de reddened) amplitudes 
are ~2-3 mJv (Jender et alllligV,: .Mirabel et al. 1998; 
Fender fc PoolevI Il998|) - 10 mJy llEikenberrv et cTl 



1998a'b' ) and ~2 mJy ((Fender fc Poo lev 200(f). Finally. 
Eiken berrv et al.l (|200flt) presented a third type of flare: 



faint (submillijansky) IR (K-band) flares whose associ- 
ation with X-ray soft-dip/soft-flare behaviour is uncertain. 

Here we present the first mid-IR (4 to 18 /im) ob- 
servations of GRS 1915+105 thanks to the sensitivity of 
ISOCAM, the IR camera on board the Infrared Space 



Observatory (ISO). These observations were part of a gen- 
eral campaign studying the interaction of the high energy 
emission from X-ray binaries with the surrounding mate- 
rial. In particular, thermal emission from heated dust was 
expected as the presen ce of dust arround G RS 1915-1-105 
had been suggested bv lMirabel et al.l (|l996(l after observ- 
in g: a reddening at th e time of a near-IR flare, and also 
bv lMartf et al.1 (|2000(l to explain the P Cygni profile of a 
He I line. In this paper we will first describe our ISOCAM 
observations and the difficulties in dereddening the data 
inherent to the mid-IR wavelength range. In section^jand 
|3|we show and interpret the results corresponding to two 
particular epochs in 1996 and 1997, when the source was 
in very different states: a possible flaring state in 1996 
and the plateau state preceeding a giant outburst with su- 
perluminal ejections in 1997 October. We also discuss the 
nature of this mid-IR emission. Section |S| summarizes our 
conclusions. 



2. Observations 

2.1. ISOCAM 

Mid-IR observations of GRS 1 915-1-105 were ca rried out 
with the ISOCAM camera ( Cesarskv et al.ll996l) on board 
the Infrared Space Observatory (ISO) satellite mission 
as part of the Guaranteed and Open Time programmes. 
Table ^ shows the log of the observations for the two 
epochs in 1996 April and 1997 October. 

2.1.1. Data reduction 

The ISOCAM data were reduced with the 
Cam Interactive Analysis software (CIA . 
iDelanav fc ISOCAM Interactive Analysis Teaml kOOd) 
version 4.0, f ollowing the standard processing outlined in 
IStarck et alJ 1x999). First a dark correction was applied, 
then a de-glitching to remove cosmic ray hits, followed 
by a transient correction to take into account memory 
effects. Pixels showing remnants of these effects were 
masked as well as side pixels insufficiently lit. The jitter 
correction was applied for 1.5" x 1.5" resolution images. 
The flat-field correction used the automatic evaluation 
except for the LWl filter where the calibration fiat field 
was used. Then individual images were combined into 
the final raster map and finally the pixel values were 
converted into milli-Jansky fiux densities. No colour 
correction was applied. 

2.1.2. Images and photometry 

GRS 1915-1-105 appears on the ISOCAM images as a 
very faint point source, at the limit of the ISOCAM 
detection for a few filters in 1996. No extended emis- 
sion was found, but the 0.3-0.6" long extensions seen by 
ISams et all (|l996b|) at 2.2 /im in 1995 July are far too 
small for ISOCAM spatial resolution. Note that these 
near-IR "jets" were suggested to be of synchrotron na- 
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Table 2. Infrared flux densities of GRS 1915+105. i^obs is the obs erved flux direc t ly me asured on ISOCAM images, 
except for the J, H and K band fluxes which were reported from iMahonev et aL I lll997l) an d tran sformed in to mJy 
units. Fdor is the resulting flux after dereddening with two different laws: from lLutz et alJ l|l996|) and from iDraind 

ET"" 
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ture and va, ri able. b ut their presence was never conflrmed. 
ISams et all l)l996a|) noted the disap pearence of the jet 
since their previous observation and lEikenberrv Sz Faziol 
lll997^ found no evidence of near-IR extended emission. 

Achieving the photometry of GRS 1915+105 was a del- 
icate task as it is located close to a very bright mid- 
IR source (J2000 coordinates a = 19'*15™11.27'' and 
6 = 10°56'27"). This so urce is the star visible in Figure 2 
of iMirabel etaP lll994l at the South- West o f the ROSAT 
error circle on this R-band image, and in iFender et al.l 
l)l997tl Figure 1 (next to star "B") where it appears as 
already quite bright in the K-band. The ^ 18" distance 
(center to center) between this star and GRS 1915+105 
corresponds to only 12, 6 and 3 pixels in our ISOCAM 
images depending on the spatial resolution. This star is 
very bright in all the ISOCAM images and its point spread 
function certainly overlap the one of GRS 1915+105. Thus 



we measured the integrated flux inside a tight circle 
around the source peak and we used an aperture correc- 
tion to obtain the observed flux i^obs- The photometry 
accuracy was estimated to ~ 30% mainly because of mea- 
surement errors due to the presence of this bright star. 



As shown in Tabled GRS 1915+105 was observed at 
different epochs with different fllters and different spatial 
resolutions which raises doubt about whether the flux vari- 
ability is real or not. Comparison with the bright source 
shows that its flux is stable on all our observations (within 
the error bars) and so the variability of the flux densi- 
ties measured for GRS 1915+105 between 1996 and 1997 
shown in Table El is real. 
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2.2. Near-IR values 

We have looked for simultaneous observations of 
GRS 1915+105 with our IS OCAM ones and we have found 
that iMahonev eraD l)l997j) observed the source in the J, 
H and Kg (2.15 /xm) bands on April 28 1996 at the same 
time as ISOCAM LWl, LW8 and LW2 filters as shown 
in Table ^ In this table the J, H and K s band fluxes 
were reported from iMahonev et al.l l)l997l) after conver- 
sion from magn itude into milliJansky units according to 
IZombec k ( 1990) reference fluxe s. They were dereddened 
using the .Lutz et alJ l)l996|) and lDraing l)l989(l laws as de- 
scribed below in section ESI These laws give at these 3 
near- IR bands very s imilar results to the more commonly 
used ICardelli et al.l l|l989t) law, so as they were used in 
the mid-IR range we kept them in the near-IR for coher- 
ence. Error bars on these ban ds were calculated bv ad ding 
observational errors given by IMahonev et al.l l^^^) ^^'^ 
errors caused by dereddening, about 10% for J and 2% for 
H and Kg, so they are respectively 12%, 5% and 4%. 

We have also looked for the maximum and minimum 
published value s in the J, H and K bands. The min- 
ima come from IChatv et al.l iim(t\. the m a xima come 
from IMahonev et all l|l997t) . iMirabel et aD l(l99fil) and 
IChatv et all l(l99fil) for the J, H and K-bands. We trans- 
formed these fluxes into milliJansky units and dereddened 
them. 



3. The dereddening problem 

3.1. Distance and visible absorption 

The distance of GRS 1915+105 is not known precisely. 
Observation of the supcrluminal e.jecta enable s us to 
constrain the maximum distance. iFender et alJ l)l999a|) 
calculated their relativistic proper motion and found 
11.2 ± 0.8 kpc as an upper limit. They also discussed the 
minimum distance, which is > 6.6 kpc. We adopt the value 
of 11.2 kpc. 

At such a distance and because GRS 1915+105 is lo- 
cated in the galactic plane with a = 19'' 15 ™ 11.5494^ and 
S = + 10° 56' 44.758" in J2000 coordinates l|Dhawan et all 
kOOd) corresponding to I = 45.3657° and b = -0.219°, 
the visible radiation from this source is heavily absorbed 
by the material along the line of sight. The optical ex- 
tinction A y was roughlv e valuat ed to be Ay ^ 30 mag- 
nitudes by iMirabel et all ll 19941) and Ay ^ 28 mag by 
iBoer et alJ l(l996|) . IChatv et all l|l996l ) derived from mil- 
limeter "'^^CO (J=l — 0) observations a total hydrogen col- 
umn density A^h = 4.7 + 0.2 x 10^^ cm~^ along the line of 
sight corresponding t o Ay = 26.5 ib 1 mag. This value was 
recently corrected bv'Chauuis fc Corbell l|2f)03[ in prep.) 
who find Ah = 3.5 ± 0.3 x 10^^ cm'^ and thus Ay = 
19.6 ±1.7 mag. In this article we adopt Ay — 20 mag. 

Note that this absorption coefficient is the main pa- 



Absorption lows 



rameter used to deredden optica l and infrared fluxe s. With 
Ay = 20 and using either the iRieke fc Lebofskvl l)l985|) 
or the Cardelli et al. ('1989') law one finds Ak = 2.2 (see 
also iChapuis fc Corbel .200,'! in prep.) and not the 3.3 



LW2 LW3 
LWl LW4 LW6 MZ-LWS LW9 




+ Mathis (1990) law 
O Dralne (1989) law 
□ Lutz et al. (1996) la 



A (^im) 

Fig. 1. Comp arison of the near-IR absorption law by 
iMathij 1I1990I) where A^ oc A "^ '' (plus sign) with lDrairi3 
l)l989l ) and lLutz et al ]fl996) laws (respectively diamonds 
and squares) in the mid-IR. The extinction coefficient is 
defined by A\/Ay where Ay is the exctinction in the 
optical V (0.55 ^m) band and Ax is the exctinction at 
wavelength A. The near-IR and ISOCAM filter wavelength 
ranges used in this article are overplotted. 



value, considered as a p ossible significant overestimate by 
iFender fc PoolevI l)l998l) but widely used in previous arti- 
cles with near-IR observations. Then for example the ob- 
served 3 mjy K-band oscillations amplitude leads when 
dereddened not to ~ 60 mJy but to ^ 20 mJy which is 
quite different from the radio oscillation amplitude. As 
a consequence the K-band flux in general is not as high 
as previously cited and thus the radio to IR spectrum 
is not as flat ijFender e t aDll997l: li^der fc Poolevll200nl) . 
altho ugh the si multa neous ra dio and millimetre observa- 
tions ( Fender fc Poolevl2000() in which there are no dered- 
dening uncertainties show a quite flat spectrum. However 
the flat spect r um a rgument may not be used against the 
Ivan der LaanI l(l966f) model of expanding clouds of plasma. 



3.2. Dereddening 

Mid-IR data need specific dereddening as the absorption 
law is highly wavelength dependent and irregular in this 
wavelength range co m pared to the n ear-IR range. The 
ICardelli et all l|l989l) or lMathisI ((l99(t laws are no longer 
valid for A > 2.5 fj,m (so in the mid-IR) as they totally 
ignore typical features such as the ones due to silicates. 
Th e Galac t ic inte rstell ar medium was studied separately 
by iDrainel l|l989l) and iLutz et all |T996) who found two 
different absorption laws, particularly in the 2.5-8.5 fj,m 
(3.5 — 12 X 10'* GH z) range. Fig. |T| shows the comparison 
betwe en the lMathis U990 ), Drai nc (1 989 ) and lLutz et"ai] 
(Il996l) laws. Accordin g to Anthony Jones (priv. comm.) 
the iLutz et alJ ((l996l ) law applies for the galactic cen- 
ter, for compact Hli regions and for particular sources 
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surrou nded by high carbon density regions. The iDraind 
l)l989|) law is supposed to apply in more diffuse regions 
so for the majority of Galactic sources. GRS 1915+105 is 
not in the Galactic center region but it is in the dense 
Sagittarius arm, and we do not know its chemical envi- 
ronment. Thus, we could not choose between one or other 
law, so we have dereddened the flux densities with both 
laws and have studied each case in the following sections. 

Note that what we call dereddened flux F^er is 
linked to the observed flux Fobs by the formula: 
^dcr ~ Fobs X 10+°'^"^^ where Ax is the extinction at 
wavelength A, calcula ted as Ay x th e ex tinction coefficient 
at A given by either iDraind 1)19891) or iLutz et all l)l996(l 
law. 

Dereddening introduces an error of about 5% into the 
flux density so it is negligible compared to the 30% error 
due to the photometry. Thus the total error on ISOCAM 
dereddened flux densities is ^ 30% of the flux densities 
(see Table I2J. 

4. A high/soft state in 1996 ? 

4.1. State of the source 

We observed GRS 1915+105 with ISOCAM for the first 
time on April 28 1996 (MJD 50201). Observed and dered- 
dened fluxes are given in Table El 

Because of the high variability of GRS 1915+105, it is 
important to know the state of this source for data in- 
terpretation. The radio and X-ray monitorings of Fig. |21 
illustrate the changes in the flux densities. The quick- look 
results provided by the ASM/RXTE team shows the in- 
tegrated count rate of GRS 1915+105 over the 2-12 keV 
band. The source behaviour was quite variable during the 
month preceeding and followi ng our IR observations (see 
also the light curve in Fig. 4 o f lPoolev fc Fendeill 997^. On 
MJD 50201, GRS 1915+105 showed a relatively high soft 
X-ray flux around 80counts/s. The corresponding BATSE 
occultation data showed no detectable hard X-ray flux in 
the 20-100 keV range on the same day. Thus, the state of 
GRS 1915+105 during our IR observations is moderately 
high and moderately soft. 

There was no simultaneous radio observation of this 
period as the Green Bank Interferometer (GBI) was 
not operating. The Ryle Telescope (RT) was observing 
GRS1915+105 at 15 GHz but not continuously as shown 
in Fig. 121 (RT plot kindly provided by Guy Pooley). There 
were observations taken on MJD 50199 and MJD 50202 
i.e. about 2 days bef ore and one dav after o ur ISOCAM 
observation (see also lPoolev fc Fendeilll99'^ . They show 
a very low flux level of about 3 ± 3 mJy classified as a 
detection (Pooley priv. comm.) so the source was in a low 
radio state through the period of our observation. 

4.2. Results 

In order to understand the nature of the IR emission from 
GRS 1915+105, we have tried to fit it using simple models. 



GRS 1915+105: 1996, 1 April - 2 June 
0.10 I 
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Fig. 2. X-ray and radio flux monitoring of GRS 1915+105 
around April 28 1996 (MJD 50201). The dotted line marks 
the begining of our ISOCAAf observation on that day. The 
upper panels show the photon flux measured by BATSE 
in the 20 — 100 kev band and the integrated count rate in 
the 2 — 12 keV band by the Rossi X-ray Timing Explorer 
All-Sky Monitor (RXTE ASM). The lower panel shows 
the corresponding flux density measured at 15 GHz with 
the Ryle Telescope (plot kindly provided by Guy Pooley). 



X-ray binaries gather many sources of IR radiation: the 
donor star, the external edge of the accretion disc, thermal 
reprocessing of the X-ray disc flux, the basis of possible 
compact jets, heated dust around the system... 

The most ob vious IR sour c e is th e donor star. It 
was identified by ICreiner et all l|2001hj) as a K-M giant 
star with a temperature estimate of ^ 4800ji5QQ ^ ^^'^ 
a magnitude of K=1 4. 5-15.0 uncorrected for extinction. 
iGreiner et all l|2001a|) estimated the Roche lobe size to be 
21 ± 4 Rq. We have fitted the star contribution with a 
black body at a distance of 11.2 kpc, with a temperature 
T and radius R. We have chosen the lowest temperature 
T = 4300 K since it corresponds to the most displaced 
spectral curve toward the mid-IR range. We have found 
the black body radius R= 19-23. 5 Rq by comparing the 
resulting black body flux in the K-band with the dered- 
dened K magnitude. Finally we have chosen R = 23.5 R© 
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Fig. 3. IR sp ectrum of GRS 1915+105 in 1996 April 28, dereddened with the (left) iLutz et alJ l)l996|) and (right) 
iDraind l)l989j) laws. The dotted lines correspond to the minima and maxima flux densities published in the J (1.25 /im 
= 2.4 X 10^ GHz), H (1.65 /xm = 1.81 x 10^ GHz) and K (2.2 /itm = 1.36 x 10^ GHz) bands. ISOCAM filters are, from 
left to right: LW3 (1.67-2.5x104 GHz), LW8 (2.5-2.8x10-* GHz), LW7 (2.8-3.53x10'' GHz), LW2 (3.53-6.0 x 10"^ 
GHz), LWl (6.0 - 7.5 x 10^ GHz) (see TableHJ. The low radio level observed with the Ryle Telescope at 15 GHz 2 
days before and one day after the ISOCAM observations was added. Tentative fits of the IR flux densities are plotted 
as indicated in the legend (see details in the text). 



in order to model the star contribution with the highest 
flux density. This black body is plotted with the dashed- 
dotted line in Fig.Elwhere it is clear that this contribution 
is faint in both near-IR and mid-IR range. 

Thus other contributions are needed to explain the IR 
emission. Near-IR fluxes can be approximatively fitted by 
a power law {Fi, (x v") with a spectral index a — 2 (thin 
solid line in Fig. but its interpretation is questionable. 
It may correspond to optically thick free-free emission but 
it would be observed at unusually high frequencies. This 
power law may also be due to the external edge of the 
accretion disc at a temperature of ~ 8000 K (and radius 
a ^ 35 Rq), it would be the low energy tail of the multi- 
black body model commonly used to account for the ma- 
jority of the soft X-ray emission in the high/soft state. 
However, as shown in Fig. |21 this power law emission is 
not sufficient to explain the mid-IR emission. Thus an ad- 
ditional component is needed in the mid-IR range. This 
component may be thermal emission from dust as plot- 
ted with the long-dashed line in Fig. 13 using a simple 
black body mod el with T = 400 K and R = 650 R^ (with 
iLutz et al.lll99^ or R = 680 Rq fwith lDrainel[il)89j) . 

The sum of the 3 components is plotted with a thick 
solid line in Fig. 13 but the result is mod erately satisfac torv 
for flux densities dereddened with the IPrainel l(l989() law 
and not satisfactory for flux densities dereddened with the 
IluIz et al.. (,1996;) law. 

4.3. Interpretation: a possible flaring event ? 

This non-satisfactory fit of the April 28 1996 (MJD 50201) 
observations may be due to a Hare occuring at the same 
moment in GRS 1915+105. 



IPoolev Fe^^ l|l 997^ described MJD 50135-50220 
as a period of varying degree of activity but with no sig- 
nificant radio emission apart from one small event. Their 
figure 4 shows indeed a quiet radio period but with a 
rather sparse coverage with the Ryle Telescope at 15 GHz. 
Nevertheless, it is possible that a brief flare or a brief 
and isolated flaring period had occurred just during our 
ISOCA M obser v ation since the J-band flux measured by 
Mah onev et alJ l)l997j) is the highest ever reported for 
this wavelength filter and since very short flares can oc- 
cur as the one seen jus t before MJD 50440 in Fig. 4 of 
IPoolev fc Fended l|l997f) . This possible brief IR Hare could 
have occurred simultaneously in the radio wavelengths but 
as it lasted less than one day, no trace is left on the MJD 
50202 Ryle Telescope observation. 

This case assume s that the near-IR observations of 
iMahonev et alJ l)l997(l must no longer be considered as si- 
multaneous. Indeed it seems unlikely when one compares 
the near-IR observati ons which lasted 12 rn inutes with 4 
minutes for each band ((Mahonev et alJl997f) with the typ- 
ical 30 minute duration of the quasi-periodic oscill ations 
of GRS 1915+105. For example. lFe"nder et alJ (|l997|) show 
oscillations with a flare duration of ~ 15 minutes, so 3 
observations of 4 minutes each taken during this kind of 
oscillation can lead to very different flux levels and con- 
sequently to a big slope, as the one between the J, H and 
Kg bands in Fig. |31 The extrapolation of this slope is not 
consistent with a detectable radio signal and this would 
be completely inconsistent with all the past recorded ra- 
dio flares corresponding to high mean flux densities of the 
order of 25 to 50 mJy. 

This possible flaring event is not identified as a pe- 
culiar variability of the flux evolution in each successive 
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image taken with ISOCAM, perhaps because of the slow 
response of the mid-IR detectors emphasized by the weak- 
ness of the source. Thus flux densities measured in the 
different ISOCAM filters during this flare have to be con- 
sidered as time averages (about 7 minutes per filter). As 
a consequence, the resulting spectral energy distribution 
does not represent a unique spectrum of GRS1915-fl05 
and one should not rely too heavily on its modelisation. 
Near-IR oscillations observed as 30 minutes dip/flare cy- 
cles in the past were considered as synchrotron emission 
(see section 1.1.). 

Another case to take into account is if no flare oc- 
curred during our ISOCAM observation on MJD 50201. 
T hen one has to explain why the J-band flux observed 
bv lMahonev et al is the highest ever recorded. A 

possible explanation would be that the accretion disc had 
been particularly bright at that time, leading to high near- 
IR flux densities coming from the outer edge of the accre- 
tion disc, but it is unlikely as the RXTE/ASM flux was 
not exceptionally high. 

Thus, the only reliable argument is that the donor star 
cannot account for most of the near and mid-IR emission 
of the source. One or several additional components 
are needed to explain this 1996 observation, such as 
synchrotron emission from a flaring event or thermal 
emission from the accretion disc external edge or from 
dust surrounding the binary system. 



5. The plateau state in 1997 

5.1. State of the source 

The second set of ISOCAM observations of GRS 1915-^105 
took place in 1997 October as noted in Table Q Fig- El 
shows the X-ray and radio flux monitoring of the source 
during this month and the begining of 1997 November. 
Because of the quasi permanent and non predictable vari- 
ability of the source in the IR, we have only plotted in 
Fig. Influx densities measured continuously on October 24 
1997 (MJD 50745). 

The GBI flux densities of GRS 1915-^105 and their 
estimated 1 a errors on this day were at about 2h30 
(UT): F2.25GHZ = 44 ± 4 mJy, Fg.acHz = 56 ± 6 mJy, 
spectral index a — 0.18 ± 0.15 {Fi, oc ly") and at about 
19h (UT): i^^2.25GHz = 44 ± 4 mJy, Fg.acHz = 40 ± 6 mJy, 
a = —0.08 ± 0.18. This last observation occurred less than 
one hour after the end of our ISOCAM observation so we 
have used it for our spectral fits (see next section) although 
we have plotted the fluxes corresponding to both hours in 
Fig. |31 On that day. X-ray fluxes show a low/hard state 
with a quite steady low RXTE/ASM 2-12 keV level of 
~ 35 counts/s and an irregular high BATSE 20-200 keV 
photon flux of ~ 0.06 ph/cni^/s. 

At that epoch, the radio and X-ray monitorings show 
that GRS 1915+105 was in a very particula r state called 
the pl ateau state (see Fig. 0] and Fi g. 1 of [ Fen der et al.' 
119994 - This state, first described bv lFoster et'al.. (,1996i) 
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Fig. 4. X-ray and radio flux monitoring of GRS 1915+105 
around October 24 1997 (MJD 50745). The upper pan- 
els show the photon flux measured by BATSE in the 20- 
100 keV band and the integrated count rate in the 2- 
12 keV band by the Rossi X-ray Timing Explorer All- Sky 
Monitor (RXTE ASM). The lower panels show the corre- 
sponding flux densities measured at 2.25 GHz and 8.3 GHz 
with the Green Bank Interferometer (GBI) and the spec- 
tral index from these two frequencies. On the latter plot 
the zero value has been marked with a solid line in order 
to emphasize the optically thick {a 0) radio emission 
of the plateau state. The dotted lines mark the begining 
of our ISOCAM observations on 1997 October. They took 
place during this noticablc plateau state preceding a giant 
radio outburst with superluminal ejections. 



and more specified by iFender et al.l l)l999a(l . results in 
an increase in the radio flux density to levels of ~ 50- 
100 mJy with an onset and a decay which can be as short 
as ~ 1 day. The plateau state is characterized by an op- 
tically thick radio emission {a ^ 0) and by a hard X-ray 
emission with a fairly high BATSE 20-100 keV flux as- 
sociated to a quasi-stable RXTE/ASM 2-12 keV flux. Its 
duration ranges from one day or two to many weeks. The 
inverted-spectrum radio emission [a > 0) corresponds to 
powerful self-absorbed quasi-continuous jets which were 
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observed bv lDhawan et alJ l)2000|l and have spectral simi- 
larities with compact jets from other b lack hole cand idate 
systems in similar hard X-ray states l|FendeJl200l|l . The 
plateau state also correspon ds in X-ray to the can onical 
low/hard C state class x of iBelloni et alJ l)2000a(l dom- 
inated by a power-law component with little or no disc 
contribution in the X-ray band. More specifically this 
term associated with the presence of compact iets ap- 
plies to xi and X3 states as studied bv iKlein-Wolt et all 
l|2002() . also called radio-lo ud or radio-plateau low/hard 
st ates bvlMuno et alJ ll200ll) and type II hard st eady states 
bv lTrudolvubovl l)200l|) . IVadawale eraD l)200l|) have even 
reported detection of synchrotron radiation in the 0.5- 
180 keV X-ray band of RXTE PGA and HEXTE during 
this state, although it has to be confirmed. 

Long (i.e. more than a few days) plateau states ap- 
pear to precede major optically thin radio flares corre- 
sponding to super luminal radio ejection events (see also 
iFender et al.ll20o3) . Plateau states are thus of further in- 
terest in order to understand the release of major superlu- 
minal outflows. As shown in Fig. 01 our ISOGAM observa- 
tion on October 24 1997 were followed by a major flare 5 
days later (on MJD 50750) with superluminal motion ob- 
served bv (Fender et al l |l999al) with MERLIN. As shown 
by the GBI spectral index, we observed GRS 1915+105 at 
the transition period from optically thick toward optically 
thin radio emission preparing the major outburst. 

5.2. Results 

The ISOGAM flux densities on October 24 1997 show 
quite different spectral ener gy distributions de pending 
on th e dereddening law used I Lutz et al]ll996l or lDraing 
Il989l) . As for the 1996 observations, we have fitted the 
mid-IR emission with very simple models given the large 
error bars on the ISOGAM flux densities (~30%). In 
Fig. Owe have plotted the contribution of the donor star 
with the same parameters: T = 4300 K, R = 23.5 R© 
(dashed-dotted line). 

As the source was in a plateau low/hard state, we have 
fitted the radio emission with a power law representing the 
synchrotron radiation from the compact jets (solid thin 
line). But whether in this state the synchrotron compo- 
nent extends to the IR range, as has been asserted by 
iFendeil 1 200lll . is questioned. This is why we have tried the 
three following cases: (1) no break in the synchrotron emis- 
sion, (2) a break occuring close to the mid-IR ran ge (re- 
spect ivel y at 9000 GHz and 31 000 GHz for the Lutz et alJ 
Il996l and lDrainelll989l laws) so that the synchrotron emis- 
sion contributes to the mid-IR flux densities, and (3) a 
break occuring before the mid-IR range so that the syn- 
chrotron emission is negligible at those frequencies. The 
optically thick part of this synchrotron component has 
been modeled using the GBI flux densities measured only 
Ih after ISOGAM thus with the corresponding spectral 
index a = —0.08. The slope of this synchrotron emission 
after the break (a = —0.8) has been fitted according to 



the ISOGAM flux densities in the LW4, LW6 and LW7 
filters dereddened with the Drainc ( 1989) law. 

Fig. El shows that ca s e (1) c an fit flux densities dered- 
dened with the IPrainel l|l989t) law as they are grouped 
around the same level (~ 20 mJy) and the a — —0.08 spec- 
tral index perfectly matches the mid-IR flux level despite 
the large error on this index (±0.18). But case (1) cannot 
fit flux densities dereddened with the iLutz et ahi l(l996^ 
law where an additional component is clearly needed 
for the LW6 (3.53 - 4.28 x 10"* GHz = 7-8.5 /zm) and 
LW4 (4.61 - 5.45 x 10"* GHz = 5.5-6.5 fim) filters. Gase 
(2) is satisfac tory for flux densities dereddened with the 
IPrainel l|l989|^ law if the break occurs in the mid-IR range 
(31 000 GHz = 9.7 //m). In all the other cases an additional 
component is needed. We then have used a black body 
model (long-dashed curves) standing for thermal emis- 
sion from heated dust, with the followi ng paramet e rs: for 
the flux densities dereddened with the iLutz et al.l lll 996l^ 
law T = 1100 K and R = 150 R© (1), T = 950 K and 
R = 330 Rq (2), T = 900 K and R = 430 R^ (3) : for 
the flux densities dereddened with the iDraind l)l989(l law 
T = 500 K and R = 850 Rg in case (3). Note that these 
values must only be considered as orders of magnitude 
since the black body model is a very simple approxima- 
tion for dust emission. 

We can summarize the results of our fits, represented 
m Fig. |S| by the sums of the components in thick solid 
lines, as the following. The synchrotron emission extended 
to the mid-IR range can account for the flux densities in 
this waveleng th range in the case of dereddening with the 
IDraind l)l989|) law. In this latter case the mid-IR flux den- 
sities can also be explained by thermal emission from pos- 
sible heated dust. Dereddening with the lLutz et al.l l(l996^ 
law implies the need for a third component accounting for 
most of the mid-IR emission, possibly thermal emission 
from hot dust. 



5.3. Interpretation / discussion 
5.3.1. Dust emission ? 

If the mid-IR emission is due to hot dust, then one has to 
explain its presence and its heating mechanism. 

First we question if the source of this dust could be 
the K-M giant donor star. Such a star has left the main 
sequence phase to evolve toward the AGB one. When 
it is isolated, its atmosphere expands and their mass 
loss increases to ~ 10^^ — 10~* Mg.yr"^ during this gi- 
ant phase. As GRS 1915+105 is a binary system, the at- 
mosphere expansion is limited by the Roche lobe size, and 
in this case the star has very likely reached its maximum 
size and has entered in the Roche lobe overflow phase with 
an accretion process characteristic of low-mass X-ray bi- 
naries. In order to produce dust, the material lost by the 
companion should escape from the black hole gravitational 
attraction, so the mass loss should exceed the mass accre- 
ti on rate. This is unlikely as the latter has been evaluated 
bv lBelloni et~l.. (,2000bil to ~ 10^* - 10~^ M© .yr" ^ during 
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Fig. 5. Rad io and ISOCAM spectrum of GRS 1915+105 in 1997 October 24, dereddened with fleffl iLutz et"al] l|l996>l 
and (right) iDrain c ( 1989) laws. The dotted hnes correspond to the minima and maxima fluxes published in the J 
(1.25 Aim = 2.4 X 10^ GHz), H (1.65 /xm = 1.81 x 10^ GHz) and K (2.2 ^m = 1.36 x 10^ GHz) bands. ISOCAM filters 
are, from left to right: LW3 (1.67 - 2.5 x 10^ GHz), LW9 (1.875 - 2.14 x 10* GHz), LW8 (2.5 - 2.8 x 10"* GHz), LW7 
(2.6 - 3.53 x 10"^ GHz), LW6 (3.53 - 4.28 x 10"* GHz) and LW4 (4.61 - 5.45 x 10^ GHz) (see TableCJ. The GBI 2.25 
and 8.3 GHz fluxes observed 14h before (higher level) and Ih (lower level) after the ISOCAM observations are also 
plotted. These latter GBI values have been used in our tentative fits of the radio to IR fluxes which are plotted as 
indicated in the legend (see details in the text). 



IR quasi-periodic oscillations even if'Bclloni et al.l ( 2^0^a^ 
has estimated it to be very low during the plateau inter- 
vals. But even if the donor mass loss was high enough, 
the wind from a giant star is supposed to be composed of 
atomic elements and not of molecular ones which are on 
the other hand needed for dust production. Moreover the 
latter is usually not enabled with isolated giant stars since 



they are too hot. Thus it is unlikely that the donor star is 
a source of dust. 

The dust could come from an external source of the 
binary system: the bright star seen in our ISOCAM im- 
ages (see section 2.1.2). This bright star is p robably in 
the AGB phase according to the study of iFelli et all 
((2000.) . In order to stay consistent with this study, we 
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took the magnitudes of this star from the ISOGAL survey 
| |http://ww w-isogal. iap.fr/): ~ 7.99 mag at 7//m (LW2) 
and ~ 7.92 mag at 15 (LW3) giving [7] ~ 7.43 mag 
and [15] ~ 7.62 mag when dereddened with Ay = 20 mag 
a nd the relat i ons A 7 = 0.028Ay and A15 — O.OlbAy as 
in lFelh et all l)2000|) . so [7]-[15] ~ 0.19 mag which corre - 
sponds to an AGB class in the figure 1 of lFelli et al.1 1|2000|) . 
If we assume that this star is located at the same distance 
from the Earth as GRS 1915+105, then the 18" angu- 
lar separation between this AGB star and GRS 1915+105 
corresponds to ~ 0.58 pc at 6.6 kpc and ~ 0.98 pc at 
11.2 kpc. An AGB star undergoes much mass loss (till 
~ 10^^ M0.yr~^) and can eject material in its surround- 
ings as far as 1 pc. But at such a distance only little ma- 
terial is provided and so it may not be able to account for 
the mid-IR emission from GRS 1915+105. 

On the other hand, if we assume that this AGB star is 
much closer than GRS 1915+105, its surrounding nebula 
caused by the mass loss probably contaminates the line 
of sight to this source. This might explain th e anom alous 
Si and Fe abundances deduced by iLee et alJ l|2002l) from 
Chandra spectra of GRS 1915+105, instead of dust that 
they suggested among other possibilities. 

The problem is also to explain the high temperatures 
(900 to 1100 K) of the dust and the distances found with 
our fits. The distance, 150 to 850 R©, represents 6 to 35 
times the chosen radius of the giant star (23.5 R©) which 
is constrained by the gravitational interactions with the 
bl ack hole (R < 21 ± 4 R© the Roche lobe size estimated 
bv lGreiner et alJl2001a|) . These parameters would not be 
su rprising for dust surr ounding some AGB stars, as shown 
bv lDanchi et al.l l)l994|) who found M stars with inner dust 
shells very close to the photospheres of the stars (3 — 5 stel- 
lar radii) and at high temperature (~ 1200 K). But for a 
gi ant star with a ma gnitude of K=14. 5-15.0 as inferred 
bv lGreiner et alJ l|2001b) for GRS 1915+105, the luminos- 
ity of the this donor star illustrated by the black body in 
dashed-dotted line in Fig. [S] is not high enough to induce 
the luminosity of the black body representing the possible 
dust on our fits. Then the dust heating in GRS 1915+105 
is unlikely to be carried out by the donor star. 

Another possible origin of the dust heating might 
be the high energy activity of GRS 1915+105. As a 
consequence, changes in the observed temperature of this 
dust between our fits in 1996 and in 1997 (when mid-IR 
flux densities is higher) might be due to changes in the 
X-ray activity. Thus the low/hard radiations coming 
from the jets or from Compton reprocessing of the soft 
disc photons by the corona would be more efficient than 
high/soft emission from the accretion disc. However, 
the efficiency of dust heatin g by X-rays is questioned 
by Ivan Paradiis et all l)l994|) who discussed the origin 
of the 10 fim emission from the low-mass X-ray b inary 
GRO J 0422+32. They found, from the study of IVoili 
lll99ll) on X-ray irradiation of interstellar grains, that this 
mid-IR emission cannot be explained by the X-ray heating 
of dust, which leads to flux densities at least two orders of 



magnitude below the observed value, unless the system is 
surrounded by a highly non-standard interstellar medium. 

Note that the hypothesis of IR emission due to dust 
surrounding GRS 19 15+105 was initially suggested by 
iMirabel et alJ l|l996ll from the observation of an IR flare 
occuring 2 days after a radio flare and corresponding 
to a reddening between the J and K band observations 
taken quasi-simultaneously. However, at that epoch the 
~ 30 minutes cycles of flares were not known and the ob- 
servations in each near-IR band took about 10 minutes. 
As a consequence these flux densities are time average 
ones and are not reliable for a spectral indication. It is 
very likely that on August 15 1 995, the observed flar e and 
reddening of GRS 1915+105 bv lMirabel eTal] l)l996|) were 
due to a ~ 30 minutes flaring event since it happened 5 
days after a giant radio outburst, now cons idered as the 
sign of super luminous ejections, and since iFender et alJ 
lll999ai I200I have observed these quasi-periodic oscilla- 
tions in radio during the decay of such giant outbursts. 

Thus if we draw the parallel with symbiotic stars (for 
a review see lMikolaiewskall2003|) GRS 1915+105 belongs 
to the ~ 80% of S-type (stellar) systems containing a nor- 
mal giant (the remaining ~ 20% are D-type containing 
Mira-type variable with a dust shell) . However it is a very 
peculiar one with a black hole instead of a white dwarf, 
with a short binary period and orbital separation leading 
to Roche lobe overflow instead of the high predominance 
of wind- accretion in usual symbiotic stars. 



5.3.2. Other: synchrotron or free-free emission 

Thermal emission from dust seems so complicated to ex- 
plain that another solution is more likely. As we observed 
GRS 1915+105 during the plateau state, which is known 
to correspond to t he pr esence of a compact jet observed 
bv lDhawan et all l)2000(l . the high-frequency extension of 
the radio synchrotron emission to the mid-I R logically en- 
sues. This hypothesis was fi rst suggested bvlFendeij (l200lh 
gathering the radio data of Poolev & Fender 7l997) with 
the only K-band flux from Bandvopadhvav et al.. (. 19*93) 
during the plateau state of 1996 August. However these 
observations were not simultaneous and Fig. |S| suggests 
that the K-band flux probably corresponds not only to 
synchrotron emission (if at all) but that the giant com- 
panion contribution may not be ncgligable in this band. 
In our study, using the Drainc ( 1989 ) law, the extension of 
the radio synchrotron emission can account for the total- 
ity of the mid-IR emission measured in several bands. It 
may be the first time that we have evidence of synchrotron 
emission in mid-IR for GRS 1915+105 in the plateau state. 
In this case, GRS 1915+105 is similar to other black hole 
candidate systems in hard X-ray stat es showing p ower- 
ful self-absorbed quasi-continuous jets llFendeill200ll) with 
synchrotron emission extending to the n ear-IR or o^ti 



cal ranges such as GX339— 4 l)Corbel fc Fe: 



■-itx or OOTi- 
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XTEJ1118+480 jMarkoff et al.ll200lHFender et alJl200l|) 
and XTE .71550-564 llCorbel et alJl200lh . 

The synchrotron emission in the mid-IR range seems 
very likely but we cannot rule out another possibil ity cor- 
respon ding to the flat spectrum obtained with the lDrain3 
l)l989(l law: optically thin free-free emission from an X- 
ray driven wind from the accretion disc. Th i s kin d of 
emission was suggested bv Ivan Paradiis et al.1 l|l994|) for 
GRO J 0422-I-32 to ex plain their observed 10 ^m flux, al- 
though iFendeil ( 200lh suggests that this detection was also 
synchrotron emission from the compact self-absorbed jet. 
The presence of a strong radia tively driven wind was also 
suggested bv lLee et al.1 1|200^ to explain ionized features 
in the X-ray spectrum of GRS 1915-1-105, although evi- 
dence of this might be difhcult to observe. Such a wind 
may be formed by X-ray heating of the disc for X-ray 
lumin osities in excess of a f ew percent of the Eddington 
limit ( Bcgcl man et al ] 11983') which is probably the case 
for GRS 1915-1-105 l(Kinel2002.) . Infrared free- free emission 
from a disc wind may be common for X-r ay binaries as 
it has already been observed in Cy gnusX-3 jPender et alJ 
ll999bHKoch-Miramond et al.l2002l) and possibly in SS 433 
llWvnn- Williams fc BecklinI Il979l iFiichs et alJ I2003L in 
prep.). 

Note that either synchrotron or free-free emission 
corresponds to t he flat spectrum resulting from dered- 
dening with the iDrain^ lll989l) law . The latter appears 



more likely than the iLutz et al.l l|l996fl law since it 
leads to satisfactory solutions to account for the mid-IR 
emission without the need for an additional component 
as thermal emission from dust (which is not easy to 
explain). Thus our 1997 observatio n enabl e s us t o prefer 
and to recommend the use of the 'Drainc fl989') law to 
deredden the mid-IR observations of GRS 1915+105. 



6. Conclusions 

We have presented for the first time mid-IR observations 
of GRS 1915+105 obtained with ISOCAM. The near-IR 
and ISOCAM observations during the high/soft state in 
1996 are difficult to interpret since they possibly took 
place during a flaring event with a too fast time variabil- 
ity compared to the response of the mid-IR detectors. The 
emission from the giant donor star alone cannot account 
for the near and mid-IR emission and one or several addi- 
tional components are needed. 

The other ISOCAM observation took place in October 
1997, while GRS 1915+105 was in the peculiar plateau 
state characterized by a flat or slightly inverted radio spec- 
trum due to the synchrotron emission from compact jets. 
This state is of further interest as it preceded a giant ra- 
dio flare corresponding to a major superluminal ejection 
event. We showed for the first time observations during 
a plateau state with quasi-simultaneous radio and mid- 
IR measurements. We discussed the different possibilities 
concerning the latter. Due to the large error-bars neither 
solution can be totally excluded but the presence of dust 



and its inferred temperature is unlikely. Synchrotron emis- 
sion can account for the flat broadband spectrum from 
the radio to the mid-IR range. Thus the mid-IR emission 
during the plateau state of GRS 1915+105 is very likely 
to be the extension of the radio synchrotron emission, al- 
though optically thin free- free emission of an X-ray driven- 
wind from the accretion disc cannot be excluded. The 1997 
observation enables us to choose between the two dered- 
dening laws: the one from IPrainel l|l989l) giving a flatter 
spectrum in the mid-IR appears more plausible than the 
iLutz et all l)l996() law. 
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